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positive masses. We investigate singular, zero temperature solutions in the Einstein 
limit but find them to be incompatible with the concept of GB gravity being a per- 
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1. Introduction 

The gauge/gravity correspondence provides a means of studying strongly coupled 
field theories via the analysis of more amenable, weakly coupled gravity theories. In 
recent years the correspondence has been applied to condensed matter systems lead- 
ing to interesting descriptions of apparent superconductors 0, |, ||. Its application 
to superconductivity is of particular interest due to the discovery of, so called, high 
temperature superconductors in the 1980 's which fall outside of the scope of current 
theories of superconductivity and are thought to be described by strongly coupled 
field theories. 

In their simplest manifestations, the gravitational duals to these superconducting 
systems consist of a black hole in an anti-de Sitter (adS) spacetime with a complex 
scalar field coupled to a U(l) gauge field. An interesting characteristic of such theories 
is that they can remain stable despite the scalar field having a negative mass squared, 
provided it is not below the Breitenlohner-Freedman (BF) bound f^. However, it is 
argued in p that the scalar field acquires negative contributions to its effective mass 
squared which, at low temperatures, can be such that the mass drops below the BF 
bound over a large enough range to render the system unstable to the formation of 
scalar hair^. This non-trivial scalar hair has a power law fall off at the boundary 
and the coefficient of this fall off can be interpreted, via the correspondence, as a 
condensate in the boundary theory. Analysis of the boundary theory has shown 

^We use the term "hair" in the spirit of the no hair theorems, |0, ^, as in something that leaves 
an imprint at radial infinity. 
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it to exhibit many of the quahtative characteristics of a superconducting system. 
There has been a great deal of recent work in this area, studying how the system 
behaves in a variety of different scenarios including altering the scalar mass and 
spacetime dimension, changing the gauge group and adding an external magnetic field 
[^, |T0|, [11], |T2|, [1^, 0, |I5], |I6|, [l^. Such systems have become known as holographic 



superconductors. 

The majority of these models are, so called, "bottom up" models^ where the grav- 
ity theory, often Einstein gravity, is, according the gauge/gravity correspondence, 
thought of as a low energy effective theory for some overarching string theory. In 
[pg, E2l pa p|, 123, Be, ETl pa pa the authors studied the stability of these systems to 



the inclusion of the Gauss-Bonnet (GB) invariant in the gravitational action which 
is believed to be the 0{a') correction to some low energy string theories [^U|, ^ 



Thus, by studying this system the authors were studying the stability of the model 
to the inclusion of higher order corrections. The findings of these papers showed 
that the qualitative features of the holographic superconductor were stable to higher 
order corrections but the details changed. These papers however, used either only 
one particular mass to focus their analysis or confined their studies to the probe, or 
non backreacting, limit. This has meant that the full panorama of the GB invariant 
on the backreacting superconductor has remained unknown. 

This paper expands upon [^, studying the fully backreacting holographic su- 
perconductor in GB gravity for a variety of scalar masses for positive and negative 
GB coupling constant a. In a recent paper, [^, it was suggested that causality con- 
straints from hydrodynamics limit the GB coupling to a G [—0.711,0.113], in this 
work, however, we permit the full range of a G (— oo,L^/4] in our study for greater 
understanding of its effect. We study the critical temperature of the superconduct- 
ing system, both analytically and numerically for a variety of scalar masses. We find 
that in a region of parameter space close to both the BF bound and the upper limit 
of a, numerical analysis is possible even at large, super-planckian scales and that its 
effect is to increase the critical temperature, not reduce it, as is its tendency in other 
regions of parameter space. The study of negative a shows to increase as a drops 
and the effect of backreaction is diminished as a gets large and negative. 

We then look at the zero temperature limit of this superconductor first showing 
analytically that there can be no regular superconducting solutions at zero tempera- 
ture for systems with a tachyonic scalar field and place strict constraints on systems 
where the mass is positive. Following |jl3| we relax the constraint of regularity for a 



system without a black hole, allowing logarithmically divergent terms. We find that 
such solutions are incompatible with the concept of GB gravity being a perturbation 
of Einstein gravity. In the absence of reliable numerical zero temperature solutions 
we then use two analytic arguments to place bounds on the critical values of the 



^For an example of some "top down" approaches see jl^, 19, pO| 
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constants of the theory about which the zero temperature phase transitions may 
occur. 

We also study the effect that the GB couphng, backreaction and the mass have 
on the conductivity of the system, finding largely quantitative not quahtative alter- 
ations, the exception being the effect of backreaction on the appearance of quasi- 
normal modes in the conductivity for a system in the vicinity of the BF bound. 
We find that increasing backreaction quickly removes the appearance of these quasi- 
normal modes within the temperature range that we are numerically able to study. 

The paper is organized as follows: Section 2 provides a brief overview of the 
holographic superconductor in GB gravity and a presentation of the equations of 
motion that are to be solved. In section 3 we study the nature of the condensate and 
in particular the critical temperature of the system for which we use both analytic 
and numerical results. In section 4 we study the zero temperature limit. In section 
5 we study the conductivity of the system. Finally we conclude in section 6. 



2. The bulk 

In this section we review the set up introduced in . This consisted of an Einstein 
Gauss-Bonnet (EGB) gravitational action coupled to a massive charged complex 
scalar field and a U(l) gauge field 



•5 = ^ / d'^xV^ 



VI a. , „abcd p A D'^b p i p2\ 

--tt + + — (^-ft riabcd — Uab + rC ) 



--F'^'Fab + |Va^ - iqAa'4j\^ - V 



(2.1) 



where = SttGs provides an explicit Planck scale in the system, g is the determinant 
of the metric, i?, Rabcd and Rab are the Ricci scalar, Riemann curvature and Ricci 
tensors respectively. The negative cosmological constant, —Q/Li^, has been written 
in terms of a length scale, L and a G {—oo, I? is the GB coupling constant. A is 
the gauge field and V' is a scalar field with charge q. 

In this paper we will study the minimal potential consisting of a single term, 
quadratic in -0 



vm 



m 



(2.2) 



where m is the mass of scalar field. 

In order to examine holographic superconductivity we look for charged, planar, 
black hole solutions^ with or without nontrivial scalar hair. We do this by using the 



^For studies of GB black holes in adS space see 33, f54l B5 
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following metric and static field Ansatze 



— {dx^ + dy^ + dz^ 



fir) LI 



(2.3) 
(2.4) 



where without loss of generality ip can be taken to be real. Lg is the effective asymp- 
totic lengthscale of this space time given by 
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The fully coupled system of gravity, gauge and scalar equations take the simple 
form^ 
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where a prime denotes a derivative with respect to r. 

In order to solve these equations we need to impose boundary conditions at 
the horizon and the adS boundary. The position of the horizon, r+, is defined by 
f{r+) = 0. Demanding regularity of the matter fields and metric at the horizon gives 



0, 



^'(r^ 



m 



/'(r+) 



(2.10) 



At the boundary we want the spacetime to asymptote adS in standard coordi- 
nates so we shall look for a solution with 



1/ — > , /(r) ~ -r^ as r — oo . 

Li 

Asymptotically the solutions of (j) and are then found to be: 



(2.11) 



(r) ~ P 
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(2.12) 



Note that there is an ij component of the EGB equations but it is not independent of (2.8) 



and (2.9) 
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where Q is the charge of the black hole (up to a factor of 47r) and A± = 2 ± 
a/4 + ■mP'Ll. We choose to set C_ = and interpret, (Ca+) = C+, where Ca+ is 
a boundary operator with the conformal dimension A_(_. If Aj. > 3 the opposite 
choice of C+ = and (Oa_) = C*- does give normalizable solutions, but will not be 
considered in this work. An example of where such a choice is made for a system 
with Einstein gravity can be found in [|rT[ . 

This paper is concerned with the effect that varying the mass has on the super- 
conductor. We shall choose a sample of masses, greater or equal to that determined 
by the BF bound, m? = — 4/Lg, in order to observe the effect on the boundary 
theory. Each choice of mass will be fixed with respect to the adS lengthscale, Le, in 
order for the dimension of the boundary operator to remain constant with respect 
to variations in a. 

In the next section we solve ( p.6[ - |2.9| ) numerically, reading this l/r^+ fall off of 
the scalar field to obtain (Oa+) for a range of temperatures given by 



(2.13) 



r=r+ 



For numerical convenience we will use the scaling symmetries of the metric, 
to set L = Q = g = 1. With this rescaling is the parameter used to vary the 
backreaction of the fields on the metric; if k"^ = 0, referred to as the probe limit, the 
fields decouple from the metric entirely. 



3. The boundary 

We wish to study the effect that the inclusion of (GB) backreaction has on the 
holographic superconductor for the full range of masses. We begin by analysing the 
operator, (Ca+), as a function of temperature, as seen in figure |I[ This plot shows 
that at high temperatures the expectation value of the boundary operator is zero 
which corresponds to the scalar field having a trivially zero profile in the bulk. As 
the temperature drops below some critical temperature, Tc, the bulk scalar obtains a 
non-zero profile which, on the boundary, is interpreted as the operator 'condensing' 
out of its vacuum. The critical temperature is a feature particular to each system 
and can simply be read off from such plots. 

Varying a, and w? produces qualitatively similar plots to those of figure |l| 
with the key differences being a variation in T^. As well as obtaining the exact value 
of Tc from such numerically produced plots, a quicker, but rougher understanding 
can be obtained from an analytically generated lower bound on Tc first introduced 
in This bound is found by studying the scalar field equation in the vicinity of 



Tc. At temperatures just below Tc, the scalar field is small, ■?/' <^ 1, and the metric 
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Figure 1: Plot of the condensate as a function of temperature; a = = 0, m? = — 3/Lg. 



and gauge field will have the form 
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up to corrections of order 0{'il)'^). Thus, in this region, the scalar field satisfies a linear 
equation, ( p.7|) , with / and cf) taking their background values. By letting Y = r^ip 
and rearranging the field equation for Y implies that if a solution exists, then the 
integral 



r+ 



fo^ LI 



3/^ 



dr = — 



/2 



dr < 



(3.3) 



is negative. For much of parameter space this integral is negative at large T, and 
positive as T — )■ 0, thus, observing where it changes sign provides a lower bound on 



Tc. As is noted in |2^, the negativity of this integral does not imply existence of a 
solution to the linearised equation near but is simply a necessary condition on one 
if it exists. Therefore, any result obtained from the lower bound must be supported 
by numerically calculated values of Tc. 

Figures § to § show both the analytic lower bound (as lines) and numerical 
values (as points) of Tc for different values of a, and m^. Figure || demonstrates 
the dependence of Tc on m^, focussing on a > 0. It is possible to find superconducting 
solutions for > 0, indeed we found solutions up to a mass of ^ 0.4 for = 0, 
reducing slightly as was increased. The findings of [p6| suggest that solutions 



exist at even larger values but that numerical solutions become difficult to obtain 
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due to an intriguing "warping" of the space of permissible boundary conditions. The 
solutions that we did obtain for a small positive mass were only marginally different 
to those of = and so these plots have not been included. In the plots that 
are shown we see that in the majority of parameter space the effect of increasing 
backreaction is to reduce the value of Tc. However, as a — )■ and — )■ — 4/Lg 
the effect of backreaction is reversed and actually increases T^. This can be explored 
in more detail by plotting as a function of k^, as seen in figure |^. This plot clearly 
shows that in this very narrow region of parameter space the effect of backreaction 
can be to increase the critical temperature of the system substantially above its 
value in the probe limit. The ability to reach super-planckian values of backreaction 
has been verified numerically up to ~ 150. It is also interesting to note that 
as one approaches this regime the lower bound on the critical temperature becomes 
significantly less accurate. 

It is straightforward to extend this analysis to a < 0, as shown in figure ^. In 
this regime the effect of altering the mass is less marked so one mass of = — 2/Lg 
has been chosen as a representative sample. These plots show that as a becomes 
more negative the critical temperature increases. Whilst this increase becomes more 
and more gradual as a is reduced it appears that an arbitrarily large Tc can be 
obtained by an appropriate choice of a. These plots also show that the effect of 
backreaction is, in all cases, to reduce Tc, but as a becomes large and negative its 
effect is diminished. This can be understood by looking at the action, (|2.1| ). In the 
Einstein limit the curvature of the spacetime scales with n. When |a| is large the 
higher order curvature terms dominate meaning the curvature scales as ^/K and thus 
the effect of backreaction on the spacetime is reduced. 

In an attempt to provide a clearer picture of the characteristics noted above we 
can use the analytically calculated lower bound on Tc and scan through the parameter 
space available to generate the lower bound on a surface of Tc as seen in figure |[ 
Whilst these plots are only a lower bound to the true surface of Tc they do exhibit 
some of the interesting characteristics of the system that have been supported by 
exact numerical results. We see immediately how altering the mass of the scalar field 
dramatically alters the nature of this superconducting system. 

4. Zero Temperature Superconductors 

There is a great deal of interest in the zero temperature limit of these supercon- 
ducting systems and in particular in the phase transitions that happen there. Most 
phase transitions are triggered by the thermal fluctuations of the system but at zero 
temperature, where there are no thermal fluctuations, phase transitions are triggered 
by the quantum fluctuations associated with Heisenberg's uncertainty principle. The 
critical points about which these zero temperature phase transitions occur are called 
quantum critical points (QCPs). It is thought that in certain regimes the effect of 
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Figure 2: Plot of Tc as a function of a for a variety of k.^; Lines represent the analytic 
lower bound, and the points represent numerically obtained values. The solid black lines 
and circular points corresponds to = 0.0, solid grey lines and square points to = 0.05, 
black (large) dashed with triangular points to = 0.2, grey (large) dashed and diamond 
points to = 1 and black (small) dashed to = b. 



the QCP can extend to finite temperature giving rise to unusual physical phenom- 
ena. For real superconducting systems it is impossible to reach the absolute zero 
temperature required to study these QCPs. However, this is not the case with our 
theoretical models leading to a great deal of recent activity in this direction. 
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Figure 3: Plot of Tc as a function of at a = 0.24999 for different masses; For each mass 
the analytic lower bounds are represented as lines and numerical values as points. The 
black solid line with circular points corresponds to rr? = —A/L'^, grey with square points 
to = — 3.75/Lg and the black dashed line with triangular points to = — 3/Lg. 
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Figure 4: Plot of Tc with a for tv? = —2/ Li. (a) shows the region a G [—1,0.25] and (b) 
shows the same plot but for a G [—100,0.25]. The lines correspond to the lower bound, 
points to numerically obtained values of Tc- The black solid lines (and circular points) 
correspond to = 0; solid grey (and square points) to = 0.05; dashed black to = 0.2 
and (smaller) dashed black (with diamond points) to ac^ = 5. 



For holographic superconductors the temperature of the boundary theory is gov- 
erned by the temperature of the black hole in the bulk spacetime. The temperature 
of a black hole in our system is given by 



T=i-/'(r)e^« 



(4.1) 



r=r+ 



In general the temperature of a black hole can approach zero in a variety of ways 
depending on the type of black hole. For example, above the critical temperature of 
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Figure 5: Plots (a), (b) and (c) show the surface of a lower bound on Tc for 
— 3/Lg and respectively. 
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our system the black holes are simply Reissner- Nordstrom black holes in GB gravity. 
Such black holes have 



4r+ 



(4.2) 



which means that the mass and charge can balance such that the temperature goes 
to zero at finite r+. This is not the case for the uncharged Schwarzschild black hole 
that arises when = 0. This has f'{r^) = Ar^/L"^ and the zero temperature limit 
is approached when r+ — )• 0. To study this limit of the holographic superconductor 
we must investigate the hairy black hole. A priori it is not immediately obvious how 
such black holes approach zero temperature; is it found at some finite r+ or when 
r+ — ?► 0. From the numerical solutions calculated above it seems that the latter may 
be true since the temperature is reduced by reducing r+ and within the range studied 



- 10 - 



there have been no apparent zero temperature solutions at finite r+. However, it is 
numerically very difficult to approach r+ = from some finite value and it is possible 
that a zero temperature finite r+ solution exists beyond the scope of the numerics. 



In [|T3] the authors calculated numerical results for holographic superconductors in 
a regime where r+ is precisely zero and the results of which, reassuringly, seemed to 
correspond the asymptote of their finite r+ solutions. However there are problems 
with this approach. Firstly, that the r+ = spacetime is not continuously connected 
to the finite r+ spacetime which introduces an element of uncertainty into the results. 
Also the results that they obtain are singular, which raises additional concerns. 

We can attempt to find information about the true nature of the zero temper- 
ature superconductor with a little investigation of the field equations ( p.6| ) to ( |2.9|) . 
In particular we will ask whether these equations permit the existence of zero tem- 
perature, regular solutions with a non-trivial scalar field? We will show that this 
is not the case. We extend the work of by showing that there are no regular 
zero temperature solutions, including those with r+ = for scalars with tachyonic 
masses. We also address scalars with rn? > 0. 

We begin by imposing that our system be regular. This will be true if the energy 
momentum tensor, T^,^, is non-singular in coordinates that are locally regular at the 
horizon, or indeed, at r = if there is no horizon. Using Eddington-Finkelstein 
coordinates defined hj v = t + r* and r = p where r* is the tortoise coordinate 
defined by 

dr* = ^, (4.3) 
the following combinations of the energy momentum tensor must be regular 

T,, = Tu = fe'T^, (4.4) 
T.p = ^ = -e^T/, (4.5) 

Tpp = Trr + ^ = j{T}-Tn. (4.6) 
^ gives the most restrictive constraint, namely that 



+ V'^ < oo (4.7) 



must be finite and hence each of the individual terms must also be finite. We wish to 
assess whether the field equations permit these constraints to hold for a non-trivial 
solution at zero temperature. The field equations are unchanged by the coordinate 
transformation and we are free to use ( |2.6|) to ( |2.9| ) in our analysis. 

Note that ( p.8|) , plus the regularity of (Tq — T^)/f implies that z^'(r+) is regular. 
If z^'(r_(.) is regular then i^{r^) is regular and e''^'^+^ ^ 0. Thus from the definition of 
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the temperature of our black hole, ( [4.1|) , the requirement of zero temperature must 
imply that f'{r^) = 0. 

We shall now study what effect this constraint has on the scalar field equation 



^2 J,2 



N" + ( - + ^' + 4 V^' + (t^ - ) ^ = 0. (4.8) 



r fj V/e'" 

The terms containing ip" and ip' go to zero at the horizon by the regularity of ip'{r^) 
and the fact that f'{rj^) = 0, thus the last term must also go to zero. This implies 
that either V'(r+) = or = If ^(r+) ^ then, by (gj), ^ ^ which 

implies that = 0. Our analysis does not rule out the existence of regular zero 
temperature solutions for this choice of mass. In fact, it seems likely that such 



solutions do exist in light of [|13|, where similar solutions were found for a system 
in four dimensional Einstein gravity. We leave the search for such solutions in this 
system to future research. 

To investigate non-zero masses we consider ip{r^) = 0. If < then all the 
leading order terms of ( |4.8| ) have the same sign and cannot balance irrespective of 
whether r+ is finite or zero. Thus there can be no regular, superconducting solutions 
at zero temperature for scalar fields with tachyonic masses. 

Turning to > 0, where ip{r+) = 0, it is possible to place strict constraints on 
these masses if solutions exist. Using the field equation for /: 



^2 / r LP' 3 



(4.9) 



we see that if 0'(r+) = then r+ = and /(r) ~ t"^ / Li as r — > 0. From ( |4.8| ) we can 
then infer that 0(0) = as otherwise q^ip'^ip / fe^^ would be the only term at leading 
order. Then from the field equation for 

+ -2g2^0 = O, (4.10) 

we see that the last term is sub-dominant and the remaining terms cannot cancel. 

If 0'(r'+) 7^ 0, ( |4.10| ) implies that r_|_ 7^ 0. Then the leading and next to leading 
order terms of (|4.9| ) give 

fi-'A (4.11) 



g2v+ ^2^2' ■>+ 12 



By using these expressions in ( [4.8|) we obtain an equation for the allowed masses at 
zero temperature 

= ^(n2 + n) + \, (4.12) 
where n > 1 is the leading power of (r — r_(_) in an expansion of ip about r = r_,_. 
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This expression shows that there can be no regular solutions for < < 24:/ L?'. 
Thus if positive mass solutions do exist they can only be found at very large m? 
and/or backreaction; substantially above the values for which finite temperature 
solutions have been found. We also see that, unlike the finite temperature system, 
the "allowed" values of rn? are directly related to k^. These observations suggest 
that this positive mass result may be spurious. 

A key result of the above analysis is that the zero temperature limit of our super- 
conducting systems with tachyonic scalars is not regular. We now wish to investigate 
this in a little more detail. In |jl3| a zero temperature solution was presented in which 
the spacetime, with no black hole, possessed logarithmic divergences as r — 0. Such 
solutions can be found for our system in the Einstein limit but it becomes clear that 
they cannot be consistent with the idea of GB gravity as a perturbative expansion of 
Einstein gravity. The reason is that the logarithmic divergences of the metric cause 
the curvature invariants, such as the Riemann tensor and Ricci scalar, to diverge at 
r = 0. Since the GB terms involve higher order combinations of these invariants 
than Einstein gravity this singular behaviour will immediately be dominated by the 
GB terms as a becomes non-zero. If this is the case the concept of GB gravity being 
a perturbative correction to Einstein gravity is destroyed and the validity of such a 
solution must be questioned. 

The manifestation of this problem on the fields themselves can be seen from a 



near horizon expansion. Following [13|, for a = 0, one can find a set of boundary 



conditions consistent with the field equations 

= 0or^(-logr)i/2 ^ e^^ = irOogr)-^ + ... (4.13) 



where /3 = — l + yl — and (po/K is free parameter that can be used to tune the 
system. This Ansatz is consistent with the field equations provided 4g^ > —rn^H? 
and after integrating the fields out from the horizon one finds the asymptotic profiles 
to be consistent with (^^.121). Unlike in the four dimensional system we were unable 



to find an appropriate value of (po/K to remove the source of the boundary operator. 
As a result these solutions do not strictly describe a holographic superconductor. 
However, they are valid solutions of ( p.6|) to (p.9|) and can be used to demonstrate 
our point. 

The problem for non-zero a arises because a appears in the equations of mo- 
tion, (|;6D to (U), like (1 - ^). From (|J) it is possible to show that if /(r) = 
/sr'*(— logr)* then s < 2 which means that //r^ has at least a logarithmic singularity 
for t > 0. Since for a = 0, t = l>0 this means turning on a immediately incor- 
porates new, singular behaviour at r = which destroys the perturbative relation 
between GB and Einstein gravity. 
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We have shown that there can be no regular solutions to our system at zero 
temperature, except possibly for massless or very massive scalars, and we have also 
given cause for caution when considering non-regular solutions. It is possible that 
consistent, non-regular, zero temperature solutions can be found that respect the 
relation between Einstein gravity and GB gravity but it seems unlikely. However, we 
can still find out information about the nature of this system in the zero temperature 
regime in the absence of such solutions. There are two analytic techniques which can 
provide bounds on the the critical values of the constants at the QCP. The first 
bound is found by simply taking the zero temperature limit of the analytic lower 
bound argument for Tc. In other words, taking the zero temperature limit of plots 
such as those in figure ^ These are shown as blue lines in figure ^ 

Another bound can be found by studying the stability of the vacuum solution 
to forming scalar hair. Following p, ^ we study a perturbation of ip about the 
vacuum solution by using the Ansatz ip = '?/'(r)e~*'^*. Assuming ip{r) ^ 1 the effects 
of backreaction can be ignored (since its effects occur at O^ip"^)) and the scalar field 
equation becomes 

with and /(r) taking their vacuum values ( ^31) and ( p.2|) . 

The system is unstable if the field equation shows this small perturbation to 
diverge. This will be the case if there is a normalizable solution to ( 4.14 ) with 
ingoing boundary conditions at the horizon such that u has a positive imaginary 
part. In general this equation can be solved numerically providing us with a bound 
on the critical values of the constants for general T. However, we are just interested 
in the T = case for which an analytic expression can be obtained. 

For zero temperature our (extremal) GB Reissner Nordstrom black hole has 

and 

1 24 
f{r) = -fl{r-r^r + ... with fl = -. (4.16) 

It is then a simple exercise to expand ( |4.14|) in the vicinity of r+. It is suggested in 
[0, ^ that since we are concerned only with the onset of an instability it is sufficient 
to consider only the threshold case of a; = 0. In this case we find the solution to the 
expanded field equation to be 

(4.17) 
(4.18) 
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If the expression inside the square root goes negative then ^p will turn imaginary 
and oscillate infinitely many times before reaching the horizon which, according to 
[Q, indicates an instability. This provides a criterion determining the onset of an 
instability at extremality. Using ( 4.15 ) and ( [4.16 ) if 



ri3 



3 + m^L^<:^^ = :^ (4.19) 



the black hole is unstable to forming scalar hair^. 

Figure |^ shows both this bound (as grey lines) and that taken from the zero 
temperature limit of the plots in figure ^ (blue lines). The regions below each of 
the curves are the regions of parameter space for which the system is unstable to 
forming scalar hair, as indicated by each bound. It was suggested in ^ that the 
grey lines are not simply a bound but actually indicate the location of the QCPs 
in the system. Assuming that the true surface of critical temperature is continuous 
these plots immediately show that this cannot be the case as in each plot the two 
bounds cross. This has been further verified by the calculation of a non-zero critical 
temperature for a system with m? = — 4/Lg, a = —4 and = 1 indicated by the 
red point in plot (a), which is outside the region of instability as indicated by bound 
( [4.19| ). The correct way to see these curves is as complimenting lower bounds on the 
critical value of function of a, being aware that the true critical values could 

be some way above these combined bounds. 

The plots do however indicate that and a do have a significant effect on the 
zero temperature limit of the system with both bounds exhibiting the opening up 
of a region of superconductivity at large as both a and approach their upper 
and lower bounds respectively. This observation fully supports that of figure ^ where 
numerically obtained values of were found at large, super-planckian backreaction. 
From figure ^ (a) we see that this is unsurprising since in this region condensation 
must occur before the temperature drops to zero. What, however, still remains 
unclear is why the critical temperature increases with backreaction here. Another 
interesting observation that can be made from these plots is that there can be no 
QCPs in the absence of backreaction. 

It is also interesting to see how these bounds relate to equation ( [4. 12]) which 

expresses the values that must take if regular, positive mass, superconducting 

solutions exist at zero temperature. Inserting ( [4.12 ) into ( 4.19|) shows that that 

systems with these masses can never be in the unstable region as indicated by bound 

( [4.19[ ) and it is only for large and negative a that they can in the unstable region 

indicated by (p.3[). This does not prove that these solutions do not exist but indicates 

that their existence may be unlikely. 

^This same bound can be found by observing that the near horizon hmit of the extremal RN black 
hole has topology 0^5*2 x S2 which has a different, less negative, BF bound to the full spacetime. 
Observing when the effective mass of the scalar field in the vicinity of the horizon is more negative 
then the adS'2 BF bound leads to precisely (4.19). 
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Figure 6: Plots (a), (b), (c) and (d) show a bound, at T = 0, on the critical value of /t^ 
as a function of a for m? = — 4/Lg, — S/Lg, — 2/Lg and respectively. The region below 
each of the lines is the region where the system is unstable to forming scalar hair. The blue 
lines were generated using the lower bound argument from ( |3.3| ) whilst the gray lines were 
generated using ( [4.19| ). The bounds continue to become less restrictive as m? increases 
above 0. The red point in (a) indicates a system with = —4/Lg, a = —4 and = 1 
for which the critical temperature was found to be Tg = 0.268. 

5. Conductivity 

The conductivity of our boundary theory, a, is calculated by studying perturbations 
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of the gauge field, A^. We set y4j(t,r, x*) = A{r)e ^^^Ci to be our perturbation and 
solve Maxwell's equation, 



A" + + -]A' + 



A = 0, (5.1) 



with the physically imposed boundary condition of only in going radiation at the 
horizon: 



A{r) ~ /(r) ^^^+ . (5.2) 



Here, T+ is the Hawking temperature defined at r = r_|_. In ||T3[ a very elegant inter- 
pretation of the holographic conductivity in four spacetime dimensions is provided. 
It involves the recasting of their gauge equation to the form of a one dimensional 
Schrodinger equation 

-A,,,+Viz)A = uj^A, (5.3) 

where ^ is a new radial parameter, a is then interpreted as a combination of the 
reflection and transmission coefficients of a wave passing through the potential V{z). 
Viewing it in this way allowed intuition from quantum mechanics to be used to 
understand many key aspects of the conductivity of their system. Unfortunately, due 
to the higher dimensionality of the system discussed in this paper, such a treatment 
has proven less straightforward. Transforming in to the form of (|5.3| ) requires 
a change of radial coordinate to dz = followed by a change of variable of A = 
r~^A. Proper treatment of this system via the Schrodinger equation requires A to 
be normalizable. Since A{r — > oo) is finite, A{r — > oo) is infinite and hence non- 
normalizable. Instead, to calculate the conductivity we shall follow and expand 
the solutions to ( |5.1| ) in the vicinity of the adS boundary, (r — oo) giving 

. a2 aoLju"^ ^ r . ^. 

A = ao + ^ + ^^log-, (5.4) 

with ao and a2 being integration constants that are used to calculate the conductivity: 



As was noted in [E^ and |]rT| there exists an arbitrariness of scale in the logarithmic 



term in (|5.4|) introduced during the holographic renormalization process |^9[. As 
a result there exists an arbitrariness of scale in the imaginary part of a. We shall 
take advantage of this fact in our numerical calculations by choosing an appropriate 
renormalization scale in order to present the characteristics of the conductivity most 
clearly. 

Plotting o"/Tc as a function of uj/T^ shows a number of key characteristics of the 
superconductor as seen in figure 0. This plot shows the real and imaginary parts of 
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Figure 7: Conductivity: A plot showing the real (solid lines) and imaginary (dashed lines) 
parts of the conductivity, a/Tc, function of u/Tc for = -2/Ll,a = 0.125, ^2 = 
at a variety of temperatures. The grey, red and blue lines correspond to temperatures of 
50%, 35% and 25% of the critical temperature respectively. The small oscillations at larger 
u are a numerical artefact. 

the conductivity for a boundary theory at = — 2/Lg, a = 0.125 and = 0. The 
first thing to note is a pole in lm{a) at w = 0. This, by the Kramers-Kronig relations 
which follow from causality, indicate the existence of a Dirac delta function in Re (a) 
at u = 0. This delta function, that cannot be picked up numerically, represents the 
infinite conductivity of the system. The plot also clearly shows the presence of a step 
in Re((T) coinciding with the global minimum of lm{a) which is interpreted as an 
energy gap in the superconductor. Following [|1^ we let the value of this minimum 
be Ug, the value of the energy /frequency gap. This plot also shows the effect that 
temperature has on the conductivity with the grey, red and blue lines corresponding 
to the measuring of a at 50%, 35% and 25% of the critical temperature respectively. 
This shows that reducing the temperature alters the plot only slightly; making the 
step and dip sharper and more pronounced, but does not change the value of Ug. 
Accessing lower temperatures has proved numerically very difficult. The absence of 
a reliable zero temperature solution means we can cast no light on what happens at 
as T ^ 0. 

As we approach the BF bound the plot behaves quite differently, as can been 
seen in figure |[ Now we see that lowering the temperature does dramatically alter 
the plot. At T = 0.5Tc the plots looks very similar to that of figure 0, but as the 
temperature drops the step and dip shift to higher u, developing distinct peaks which 
turn into poles. These poles are interpreted as quasi-normal frequencies that have 
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Figure 8: Conductivity: Plots showing the real (left) and imaginary (right) parts of the 
conductivity, a, as a function of lo/Tc for = — 4/Lg, a = 0.125, = at a variety of 
temperatures. The grey, red and blue lines correspond to temperatures of 50%, 35% and 
25% of the critical temperature respectively. 

moved to the real axis from elsewhere in the complex plane |Tl|, 0, As the 
temperature drops further more poles appear at higher values of u/Tc (not shown). 
It is suggested in ^T|] that in the probe limit of Einstein gravity, as T — i- the 



number of these poles diverges. Since such low temperature analysis is outside the 
scope of this paper, we can shed no light on whether or not this occurs away from 
the Einstein limit. 

We are interested in observing the effect that varying a, and has on these 
phenomena. We will begin by looking at the first case; away from the BP bound 



where temperature dependent effects are less prominent. In p2| the authors studied 
the effect of a, for = — 3/Lg. They found that increasing a above the Einstein 
limit increased the value of Ug and made the step and dip more pronounced. The 
effect of increasing was to smooth out the features of the plot but not affecting the 
value of Ug, that is until the smoothing removes the presence of the hard gap^, at least 
within the temperature range studied. Studying the conductivity for larger masses 
we see very similar results with quantitative differences rather than qualitative. The 
key information has been captured on a plot of Ug against Tc as seen in figure ^ 

The grey points in left plot in figure ^ correspond to the probe, Einstein limit of 
the superconductor. One can see that for the range of masses presented, the points 
all fall close to the line Ug = 8Tc. This observation contributed ammunition to the 



speculation, |jTT|, that this may be a universal relation. This plot shows that such a 



relation is unstable to higher curvature corrections as found in [p2| . 

The plot shows that increasing a increases Ug and largely reduces Tc, except for 



^i.e Re((T) no longer is zero for small lu. 
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Figure 9: The left plot shows against for = and vr? = 0, black (triangular) 
points; rr? = — 2/Lg, red (square) points andm^ = — 3/Lg, blue (circular) points. The right 
plot shows iOg against Tc for = — 2/Lg for = red (square) points and = 0.05 green 
(circular) points. In both plots from top to bottom the points correspond to a = 0.24999, 
0.1875, 0.125, 0.0625, 0, —0.25, —1, —10, with the grey points corresponding to Einstein 
gravity. The dashed lines have been added to guide the eye. The straight line corresponds 

to Ug = 8Tc. 



very close to a = 1/4. This has the effect of moving the point decidedly off the line. 
Decreasing the mass from = increases Ug and Tc with the greatest differences 
occurring towards the upper bound of a where variations in Tc are more pronounced. 
The right hand plot shows the effect of backreaction. Increasing has very little 
effect on Ug with the majority of the effect coming from the reduction in T^. As a 
gets large and negative the points converge corresponding to the diminished effect 
of backreaction in this regime that was noted above. We were unable to extend 
these curves to much larger negative coupling as the numerical artefacts that arise 
in our calculation of the conductivity began to obscure the key features of the plots. 
However, since the calculation of the condensate seems possible for arbitrarily large 
negative a, one might expect these curves to continue towards the axis without ever 
reaching it. 

We now turn our attention to systems at the BF bound, and in particular what 
effect a and have on the development of the quasi-normal modes. Figure 10 shows 
Re(o"), measured at T = Tc/4, for = — 4/Lg, k,^ = for a variety of values of 
a. This plot shows that increasing or decreasing a does not seem to hinder the 
development of these quasi-normal modes. The dominant effect of, say, increasing 
the the GB coupling constant is to shift the poles to higher u/Tc. This increase with 
a is particularly marked as you approach the upper limit of the coupling constant. 

Figure 11 shows the effect that backreaction has on the development of the 
quasi-normal modes with a plot of Re(cr), measured at T = Tc/4, for = — 4/Lg, 
a = 0.24999 for a variety of k^. We see that turning on backreaction very quickly 
removes the appearance of the poles, at least at this temperature; it is still quite 
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Figure 10: Plot showing Re((T) measured at T = Tc/4, = — 4/Lg and = for a 
range of values of a. From left to right: red, a = —100; blue, a = —1.0; green, a = —0.25; 
grey, a = 0; purple, a = 0.125 and black, a = 0.24999. The small oscillations are a 
numerical artefact. 
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conceivable that they may appear as the temperature is dropped. Analysis of this 
phenomenon at much lower values of a show the existence of quasi-normal modes up 



to much higher values of /t^, supporting the observation that the effect of backreaction 
diminishes as a is reduced. 



6. Conclusion 

The aim of this paper was to explore the dependence of the fully backreacting Gauss- 
Bonnet holographic superconductor on the mass of the scalar field. We began by 
studying the critical temperature, Tc, of the system. We found that in the majority 
of parameter space the effect of backreaction is to reduce T^. but that in a narrow 
region where m? — )■ —A/L^ and a — ?■ its effect is reversed and actually increases 
Tc. In this regime large, super-planckian values of backreaction are numerically 
attainable. We also found that as a becomes large and negative Tc increases and 
the effect of backreaction is diminished as the gravitational action is dominated by 
the higher curvature terms. Again, this provides a regime where large values of 
backreaction are attainable. 

We studied the zero temperature limit, proving that the system does not permit 
regular solutions with a tachyonic scalar field and place strict constraints on positive 



mass solutions. Following [|T3[ we relaxed the regularity constraint for a system 
without a black hole and permitted the fields to have logarithmic divergences. Such 
a solution was found in the Einstein limit but was shown to be inconsistent with the 
idea of Gauss-Bonnet gravity as a perturbative expansion of Einstein gravity. These 
findings show that a satisfactory zero temperature solution to this system has not 
yet been found and raises questions as to whether one can be found, except possibly 
if the scalar is massless. 

We also studied the conductivity of the system. In the region away from the 
BF bound we saw that the prominent effect of, say, increasing a was to increase 
the frequency gap, Ug. In the vicinity of the BF bound, the effect of a was to shift 
the location of the quasi-normal modes that appear there, again shifting them to 
larger a; as a was increased. Otherwise a did not seem to effect their development. 
The effect of backreaction was more notable: increasing backreaction away from the 
probe limit quickly prevented the appearance of these quasi-normal modes, at least 
within the temperature range that we were able to study. 

The findings of this research suggest a number of interesting avenues for future 
research, most notably in relation to the zero temperature limit. The first thing to 
do might be to either find, or disprove the existence of, positive mass solutions in 
this limit. If they can be found it would be very interesting to study the quantum 
phase transitions that may happen there. 

Having shown that this system does not permit regular, superconducting, zero 
temperature, tachyonic solutions, it would be interesting to see if there exist non- 
regular solutions that are compatible with the perturbative relation between Einstein 
and GB gravity. If so, it would be interesting to find out how this non-regularity can 



-22- 



be interpreted. It would also be interesting to study the conductivity of this system 
in such a case, particularly in light of which suggests that an infinite tower of 
quasi-normal modes may appear on the real axis as the temperature goes to zero. 

Another interesting line of enquiry would be to further investigate the regime of 
positive scalar masses at finite temperature. In this work we were only able to find 
solutions for very small which proved to be only marginally different to those 
of = 0. In [Q, which studied the > case in the probe limit, the authors 
found solutions at larger masses but found that solutions become difficult to attain, 
possibly due to an observed warping of the solution space which was dramatically 
enhanced by increasing m^. It would be interesting to see how the inclusion of higher 
curvature terms and backreaction may affect this phenomenon. 

Another interesting outcome of this paper is the observation that there exists two 
regions of parameter space where the numerical system is substantially more stable 
to the inclusion of backreaction. Since the inclusion of backreaction is an important 
and often complicated aspect of these systems, these regimes may provide testing 
grounds where the numerical analysis of backreaction is somewhat easier. 
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